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ABSTRACT 

Richardson,  M.D.  and  Young,  D.K.,  1980.  Geoacoustic  models  and  bioturbation.  Mar. 

Geol.,  38:  205-218. 

Two  types  of  geoacoustic  models  are  used  to  describe  relationships  between  physical 
and  acoustic  protterties  of  unconsolidated  marine  sediments;  simple  predictive,  and 
theoretical  models.  Simple  predictive  models  use  the  apparent  correlation  between  physical 
properties  of  sediments  (usually  porosity,  grain  size  or  density)  and  acoustic  properties  of 
sediments  to  predict  velocity  and  attenuation  of  compressional  and  shear  waves,  Theoretical 
models  use  elastic  properties  of  sediments  (rigidity,  compressibility,  etc.)  to  calculate 
acoustic  properties. 

It  has  been  shown  that  bioturbation  by  benthic  animals  profoundly  affects  the  physical 
properties  of  marine  sediments.  Activities  including  burrowing,  ingestion/digestion/defeca- 
lion,  tube  building,  biodeposilion,  cementation  and  metabolic  activities  modify  porosity, 
grain  size,  density,  fabric,  rigidity  and  compressibility  of  sediments. 

We  hypothesize  that  bioturbation  by  benthic  animals  alters  acoustic  properties  of 
unconsolidated  marine  sediments.  An  example  is  provided  for  predicted  effects  of  bio¬ 
turbation  on  selected  acoustic  and  elastic  properties  of  a  silty  mud  sediment.  Knowledge 
of  effects  of  bioturbation  on  unconsolidated  marine  sediments  may  increase  our  under¬ 
standing  of  the  relationship  between  predicted  and  measured  acoustic  properties, 

INTRODUCTION 


Two  types  of  geoacoustic  models  are  used  to  describe  relationships 
between  physical  and  acoustic  properties  of  unconsolidated  marine  sediments: 
simple  predictive,  and  theoretical  models.  Simple  predictive  models  use  the 
apparent  correlation  between  physical  properties  of  sediments  (usually 
porosity,  grain  size  or  bulk  density)  and  acoustic  properties  of  sediments  to 
predict  velocity  and  attenuation  of  compressional  and  shear  waves.  The 
correlations  are  usually  presented  in  graphic  form  or  by  regression  equations. 
Porosity,  density  and  grain  size  are  also  used  to  predict  bottom  loss  of  com¬ 
pressional  waves  at  normal  incidence  to  the  sediment— water  interface. 

Although  disagreement  exists  as  to  which  elastic  or  vi.scoelastic  model 
should  be  applied  to  marine  sediments  (Hamilton,  1971a,  1974a;  Stoll,  1974, 
1977),  there  is  no  disagreement  that  physical,  acoustic  and  elastic  properties 
of  marine  sediments  are  theoretically  related  and,  once  measured,  can  be  used 
to  compute  one  another. 
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Any  process  which  controls  or  tilters  any  of  the  physical  or  elastic  prop¬ 
erties  also  changes  the  acoustic  properties  of  sediment.  Bioturbation  (bio¬ 
logical  modification  of  sediments)  is  such  a  process,  and  therefore  may  be 
important  theoretically  and  heuristically  in  geoacoustic  models. 

It  is  the  purpose  of  this  paper  to  review  relevant  empirical  and  theoreiical 
relationships  between  physical  and  acoustic  properties  of  marine  sediments, 
as  well  as  the  effects  of  bioturbation  on  physical  properties  of  sediments,  in 
order  to  develop  a  working  hypothesis  on  effects  of  bioturbation  on  elastic 
and  acoustic  properties  of  marine  sediments. 

RELATIONSHIP  BETWEEN  Af’OnSTlC  AND  PHYSICAL  PROPERTIES  OF  SEDLMENTS 

Prediction  of  acoustic  properties  of  marine  .sediments  from  measured, 
extrajiolated  or  predicted  physical  properties  of  sediments  is  central  to  the 
development  of  geoacoustic  models  (Hamilton,  1971b,  1974b).  Empirical 
relationships  among  physical  and  acoustic  properties  of  sediments  have  been 
develo>'>ed  from  concurrent  measurements  of  those  properties  under  labora¬ 
tory  conditions,  in  situ  by  diver-held  or  submersible-held  probes  or  from 
freshly  collected  piston,  gravity  or  box  cores. 

The  first  empirical  equations  were  developed  in  the  1950’s  between  com- 
pressional  wave  velocity  and  sediment  median  (or  mean)  grain  size,  porosity 
and  bulk  density  (Hamilton,  1956,  1963;  Hamilton  et  al.,  1956;  Sutton  et 
al.,  1957;  Nafe  and  Drake,  1957,  1963;  and  others).  The  most  recently  pub¬ 
lished  regression  equations  between  compressional  wave  velocity  and  the  above 
mentioned  physical  properties  are  those  of  Hamilton  (in:  Morris  et  al.,  1978). 
Different  regressions  were  presented  for  sediments  from  the  continental 
terrace  (shelf  and  slope),  abys.sal  hill  (pelagic)  and  abyssal  plain  (turbidite). 

In  carbonate  sediments,  porosity  and  density  should  not  be  used  to  predict 
compressional  wave  velocity  (Morton,  1975)  and  equations  for  mean  grain 
size  versus  compressional  wave  velocity  for  continental  terrace  sediments 
apply  (Johnson  et  al.,  1977).  Attenuation  of  comiiressional  waves  in  surface 
sediments  is  estimated  from  the  equation; 

Q  =  kf"  ( 1 ) 

where  the  exponent  n  is  1,  a  is  attenuation  in  dB/m,  f  is  frequency  in  kHz 
and  k  is  constant.  The  constant  k  is  determined  for  different  porosity  and 
mean  grain  size  of  sediments  from  regression  equations  given  by  Hamilton 
(1972,  1976a). 

Prediction  of  shear -wave  velocity  in  surficial  sediments,  given  porosity  or 
mean  grain  size,  is  more  difficult  and  tenuous  because  of  the  relatively  few 
in  situ  concurrent  measurements  of  shear-wave  velocity  and  sediment 
[ihysical  properties  (Hamilton,  1976b).  The  development  of  an  in  situ  .shear- 
wave  velocimeter  by  Shirley  and  Hampton  (1978)  should  increase  our  i)re- 
dictive  capability.  Hamilton  (1971a,  1971b)  predicted  shear-wave  velocity 
in  surficial  sediments  by  first  constructing  bulk  density  versus  .sediment 
bulk  modulus  (K)  and  porosity  versus  sediment  bulk  modulus  graphs.  Recent 
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regression  equations  for  these  relationships  are  given  in  Morris  et  al.  1 197h) 
for  continental  terrai'O  sediments  and  abyssal-hill  and  abyssal-plain  sediments 
combined.  Predicted  sediment  bulk  modulus  together  with  measured  or 
[iredicted  sediment  bulk  density  (p)  and  compressional  wave  velocity  I  ) 
are  used  to  predict  rigidity  or  shear  modulus; 

^=3/-\(pV,;-  -k)  (2) 

Velocity  of  shear  waves  is  then  calculated  from  bulk  density  and  rigidity; 

=  Ci) 


Sediment  bulk  modulus  can  also  be  predicted  if  sediment  bulk  density  and 
compressional  wave  velocity  are  known  (.Morris  et  al.,  1978).  Prediction  of 
shear-wave  attenuation  from  jiorosity,  mean-size  or  bulk-density  data  is  not 
jiossible  with  reasonable  accuracy  at  this  time  (Hamilton,  1976c). 

Rayleigh  reflection  coefficients  and  bottom  loss  of  compressional  waves 
at  normal  incidence  to  the  sediment— water  interface  can  be  calculated  given 
impedance  of  both  sediment  and  water  (Hamilton,  1970).  Imiiedance  of  the 
sediment  can  be  calculated  from  bulk  density  and  compressional  wave 


velocity  or  (iredicted  by  the  empirical  relationship  between  porosity  and 
impedance  or  bulk  density  and  imiiedance.  Hamilton  ( 1970)  presents  thi'se 
empirical  regressions  for  continental  terrace,  abyssal-hill  and  abyssal-iilain 
sediments.  Rayleigh  reflection  coefficients  (/?)  me  then  calculated  from  tln> 
impedance  mismatch; 


where  V,  is  the  impedance  of  the  water  and  pjV;  is  the  impedanci'  of  tlu' 
sediment;  bottom  loss  iBL)  is  calculated  from  the  Rayleigh  n'fh'ction  co¬ 
efficient; 


BL  =  -20  log  R 


Rayleigh  reflection  coefficients  and  bottom  loss  can  also  be  predicted 
directly  from  jiorosity  or  density  values  using  regn'ssion  t'quations  given  liy 
Hamilton  ( 1970), 

Physical  and  acoustic  {irojicrties  of  unconsolidated  marine  sc'dimi'iits  are 
theoretically  related  by  various  elastic  or  viscoi'lastic  models  (Hamilton. 
1971a,  197‘la,  b,  Stoll,  1974,  1977).  This  section  reviews  some  of  th(>  theo¬ 
retical  relationships  between  physical,  elastic  and  acoustic  properties  of 
sediments  to  betU'r  understand  effects  of  bioturbation  on  those  propm'ties. 
For  the  purposes  of  this  pajicr,  Hamilton’s  ( 1971a,  1972,  1974b)  viscoelastic 
model  is  used.  Hamilton  (1971a)  show(>d  when  energy  damping  is  low,  that 
one  is  justified  in  using  Hookean  elastic  equations  to  compute  various  elastic 
constants  (e.xcept  for  attimuation  wlu're  the  viscoelastic  etiuation  must  be 
used).  Pertinent  mathematical  relationships  include  eqs.  1 .  2.  6,  6  and  7; 


r. 


(A'  t  1/3  u )'  * 
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(6)  Code* 
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A' =  p(V;-  -4/3  (7i 

\vli('r('  I’j  is  the  vehicity  of  conipn'ssional  waves  (ni  see).  A  is  the  hulk  inoiiiilu 
I  10"'  dynes/em- ),  n  is  the  shear  modulus  ( 10"'  d.\  ties/em’ ).  p  is  the  sedmient 
density  (g/em  ’)  and  V’.,  is  the  v(‘loeity  of  shear  waves  (m/sec;). 

In  uneonsolidatf'd  marine  sediments,  shear  modulus  (dynamie  rit;idity). 
shear  stnm^'th  and  attenuation  are  all  rtdated  to  tlie  same  sediment  fartors: 
friction  in  sandy  sedimimts  ami  cohesion  in  silts  and  clays  (Hamilton.  1  071a). 
In  sand,  the  friction  hi'tweeii  t,'rains  is  a  function  of  the  numher  of  interf^rain 
contacts  (related  tc)  porosity,  imsm  grain  size,  sorting,  compiict ion ).  the  sur* 
face  area  of  grains,  and  the  interlocking  of  grains  (related  to  angul.-irity  i. 
iiigidity.  siieai  strength  and  attenuation  are  increased  by  lower  (lorosity 
(givi'ii  constant  graiti  size),  decreased  grain  siz('.  increased  packtng  and  greater 
angularity  ( Hamiltcm,  1972).  In  silts  ;md  clays,  cohesion  is  relati'd  to  tlie 
distance  lietween  partu  les,  inti  rpartiele  bonding  and  attraction,  sediment 
faliric,  clay  miiK'ralogy,  concentrtitions  of  inti'rstitial  I'li'ctrolytes.  organic 
composition  ;md  ttepo,sitional  history  (Hamilton.  197  1;k  Young  and  Southard. 
1978).  Iiigidity.  shear  .strength  and  attenuation  art*  lower  in  high-porosity. 
>mall-grain  sized  silts  and  clays  due  to  the*  reduction  in  strength  of  interjiart iek 
attractive  force's  and  ri'duced  interptirticle  contiiet.-'  (Hamilton.  1972), 

In  spite  of  the  obvious  relationships  betwi'cn  shear  modulus,  shear  strength 
and  atti'iiuation.  they  are  not  used  to  cominite  oni'  another.  Dynamic  rigidit\’ 
is  computed  from  density  ;ind  she;ir-\vave  vehicity.  shear  strength  is  measured 
by  viseometer.s  or  peiK'trometers.  and  attenuation  is  estimati'd  from  porosit\' 
and  mean  grain  size.  Sliear-wave  vt'locity  is  directly  rekited  to  rigidity  and 
density.  .\ny  reduction  in  rigidity  results  in  a  reduction  in  shear-wa\'('  velocity. 
Shear-wave  velocity  is  also  ri'lated  to  the  sanu*  frictional  and  cohesive  factors 
controlling  shear  modulus,  shi'ar  stri'iigth  and  at  tenuation. 

The  Inilk  modulus  or  incompressibility  of  sediments  is  com|)osed  of  three 
separate  bulk  moduli:  ( 1 )  the  bulk  modulus  of  pore  water.  ( 2)  the  Inilk 
modulus  of  mineral  grams,  and  (3)  bulk  modulus  of  thi*  frame  (Hamilton. 

1  971a  I  and  sediment  porosity.  The  bulk  modulus  of  pore  watc'r  and  of 
mineral  grains  varies  litth*  in  marine  sediments,  therefore,  thi'  contrilnit ion 
of  pore-water  modulus  and  minertil-grain  modulus  to  the  seiliment  bulk 
modulus  is  primarily  dependent  on  porosity.  The  frame  bulk  modulus  is 
dependent  on  properties  and  factors  influencing  sediment  fabric,  e.g..  porosity 
grain  ^ize,  pure-wati'r  composition  and  depositional  history.  The  si'diment 
bulk  modulus  increases  with  increasing  density,  decreasing  porositv’  and 
increasing  grain  size  (  Hamilton.  1971a). 

Compressional  wave  velocity  is  dependent  on  mcompresNibilitv.  rigtdilv'. 
and  density  of  the  seflimenl.  Since  the  rigidity  of  mo>l  unconsolidated 
marine  sediments  is  low  and  si-diment  ilensity  varies  little,  changes  m  mcom 
pressibility  tend  to  dominati'  compressional  w.ive  velocity 
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EFFECTS  OF  BIOTURBATION  ON  PHYSICAL  PROPERTIES  OF  SEDIMENTS 

Measurements  of  physical  properties  of  unconsolidated  marine  sediments, 
such  as  those  for  porosity,  bulk  density,  mean  (or  median)  grain  size,  co¬ 
hesion  and  compaction,  are  indirect  indicators  of  jirimary  properties  of 
mineral  composition,  sizes,  shapes,  orientations  and  packing  of  com))onent 
grains,  as  well  as  characteristics  of  gaseous  and  fluid  substances  filling  inter¬ 
stitial  spaces  among  the  grains.  Indirect  geotechnical  measures  of  primary 
properties  of  sediments  are  useful  for  practical  engineering  applications 
(Valent,  1974)  and  have  heuristic  value  regarding  behavior  of  sediments 
upon  being  perturbed  (Keller,  1974). 

It  is  well  known  that  marine  organi.sms  drastically  affect  all  three  phases 
of  sediments:  the  gaseous,  liquid  and  solid.  For  exami)le,  methanogenic 
bacteria  can  produce  sufficient  methane  to  saturate  pore  water  and  form  gas 
bubbles  resulting  in  an  acoustically  impenetrable  bottom  (Schubel  and 
Scheimer,  1973).  Further,  numerous  microorganisms  (Fenchel,  1969)  and 
deposit-feeding  invertebrates  (.Aller,  1978)  have  been  shown  to  produce 
radical  changes  in  pore-water  chemistry.  Marine  organisms,  including  bivalve 
molluscs,  crustaceans,  and  fish,  are  known  to  alter  clay  mineralogy  by  their 
digestive  processes  (Anderson  et  al.,  1958;  Pryor,  1975).  It  is  beyond  the 
scope  of  this  paper  to  address  the  gaseous,  liquid  and  solid  phases  directly  as 
they  are  affected  by  marine  organisms.  Therefore,  we  are  restricting  this 
discussion  to  affects  of  benthic  (bottom-dwelling)  animals  upon  the  bulk 
properties  of  sediments,  specifically  related  to  sediment  fabric  as  influenc¬ 
ing  acoustic  properties. 

Sediment  fabric  can  be  defined  simply  as  the  spatial  arrangement  of  com¬ 
ponent  sedimentary  particles.  The  relationship  of  grain-to-grain  contacts  is 
partially  understood  for  silicate  and  carbonate  sands  but  is  poorly  under¬ 
stood  for  cohesive,  fine-grained  sediments  (Richards  et  al.,  1976).  Even  a 
fundamental  parameter  of  such  great  significance  as  the  in  situ  size  of  fine¬ 
grained  sediments  is  ))oorly  known.  In  those  in.stances  where  attempts  have 
been  made  to  measure  water-stable  aggregates  of  silty  muds,  rather  than 
completely  disaggregating  particles  prior  to  analysis  (Folk,  1974),  it  has  been 
shown  that  a  major  portion  of  the  clay-  and  silt-sized  particles  are  bound  into 
sand-sized,  low  bulk-density  fecal  pellets  (e.g.  Moore,  1931;  Rhoads,  1967; 
Young,  1971).  Sediment  fabric  of  clays  is  influenced  by  a  number  of  factors, 
including  clay-mineral  composition,  ionic  composition  of  pore  water,  amount 
of  pore  water,  concentration  of  organic  matter  and  overburden  pressure 
(Meade,  1966).  Few  studies  in  which  fabric  of  relatively  undisturbed,  uncon¬ 
solidated,  marine  sediments  from  nature  have  been  observed  via  electron 
microscopy  indicate  a  “loose,  open,  random  arrangement  of  particles” 
(Bowles,  1968a).  The  influence  of  benthic  animals  upon  preferred  orientations 
of  sedimentary  particles  is  not  well  understood,  nor  are  the  altered  states  of 
sediments  resulting  from  bioturbation  fully  comprehended. 

The  word  “bioturbation”  (Richter,  1952,  in;  Schafer,  1972)  from  pale¬ 
ontology  is  used  in  this  paper  in  the  all-inclusive  .sense  of  Rhoads  (1967). 
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i.o.  “biogenic  reworking”,  as  “all  biologic  activity  physically  altering  sedi¬ 
ment”.  Bioturbation  has  been  shown  to  influence  the  following  properties 
of  sediments,  among  others:  (a)  porosity,  median  grain  size,  and  bulk  density, 
(b)  compaction  and  cohesion,  (c)  particle  orientation  and  distribution,  and 
(tl)  microtopography.  I'hese  properties  are  discussed  separately. 

Porosity,  median  grain  size  and  bulk  density 

Porosity  and  median  grain  size  increase  markedly  and  bulk  density  decreases 
as  a  result  of  bioturbation  (Rhoads  and  Young,  1970;  Young,  1971).  Increa.sed 
(lorosity  and  median  grain  size  result  (irimarily  from  increasing  overall  pore 
space  and  sediment  [larticle  size  by  th('  binding  of  small  constituent  [larticles 
of  clay-  and  silt-sized  si'diment  into  sand-sized  fecal  (xdlets  by  benthic 
animals.  Bulk  density  decreases  as  a  result  of  greater  proportion  of  low- 
density  pore  water  to  higher-density  sediment.  The  fecal  (M'llets  themselves 
may  have  high  porosity',  i.e.,  (Risk  and  Moffat.  1977).  Highly  Ihotur- 
bated  muds  may  reach  ‘d'ZVr  ftorosity  (Rhoads,  197  1).  The  u|)j)(‘r  several 
centimeters  of  aby.ssal  sediments,  which  are  intensively  reworked  by  deep-sea 
benthic  animals,  range  from  60  to  907;  porosity  (Hollister  et  al.,  1975). 

Compaction  and  cohesion 

Comiiai'tion  and  cohesion  of  marine  sedimimts  ari'  measured  in  situ  Iry 
rotational  viscometers  (e.g.  T-bar  spindles,  vane  shears)  and  Iry  pcandrometers. 
These  measures  are  directly  relatc'd  to  one  another  and  indirectly  related  to 
(lorosity  and  l)ulk  density.  Cohesive  muds  hav('  high  erosioiial  resistance 
because  of  electrochemical  interi)article  bonding  and  organic  binding  (\’an 
01|)he;i,  1966).  Highly  bioturbated,  i.e.,  [relleuzed.  muds  at  the  sediment- 
water  interface  have  little  cohesive  strength  (Bokuniewicz  et  al..  1975)  and 
are  easily  eroded  (Rhoads  and  Young,  1970).  Shear  strength  of  a  silty  muii 
sediment  has  been  shown  to  incr<‘ase  from  10  g/cm'  to  more  than  IS  g enr 
as  successive  vane  shear  measuremi'iits  w(>re  taken  20  cm  away  to  immeiliati  ly 
adjacent  to  a  membranous  tube  of  a  burrowing  aiUMnone  (  Rowiy  197  1). 

Myers  ( 1977)  has  demonstrat('d  both  increases  and  decreases  in  shear  strength 
of  a  sandy  s(>diment  resulting  from  bioturbation  by  two  different  spei  ies  of 
benthic  animals.  Rii'hardsand  Parks  (1976)  have  altributi'd  high  variance  of 
shear-strength  measurements  of  sediments  in  a  Norwi'gian  fjord  to  effects  of 
bioturbation.  Clearly,  the  resulting  fabric  of  sc'diments  at  any  given  time 
n>[)re.sents  a  balance  of  predominating  influences  of  bioturbation  by  liiffereiit 
species.  .An  open  fabric  from  ingestion/defei-ation /pellidizing  activities  of 
mobile  deposit  feeders  results  in  low  eomiiaction  and  cohesion,  whereas  a 
com|)acted,  cohesive  and  tight  fabric  results  from  extensive  burrowing  and 
tube-building  (  Rhoads,  1971 ). 


'  .All  s<'<linifnt  Wilier  Cdnlenl  vulues  .ire  presenter!  as  perei'nt  porusily  in  this  paper  VVliere 
necessary,  percent  water  content  (in  sensu  Hhoarls,  1971.  or  Keller,  1971)  were  recalculated 
■issuminn  ii  mmeriil  uriiin  density  of  2. tint)  u/cm'  and  a  water  ilenslly  of  1.021  g  cm  ' 
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Particle  orientation  and  distribution 

BioUirbation  prodiK'(>s  oither  homogonoous  fabric,  i.c.,  random  particle 
orientation  and  random  size  distribution,  or  heterogeneous  fabric,  i.e., 
aggregations  of  grain  sizes  and  compositions  (Rhoads,  1971).  As  previously 
discussed,  effective  i)article  size  may  be  increased  by  pelletizing  activities.  In 
addition,  individual  grains  may  tx'  .selectively  ordered  and  segregated  within 
fecal  pellets  in  terms  of  size  and  composition  (Pryor.  1975).  Selective  feeding 
upon  .small  particles  by  de[)osit-feeding  animals  results  in  graded  bedding 
(Rhoads,  1967),  X-radiogra))hic  opaque  layers  of  denser  grains  (Rhoads  and 
Young,  1970)  and  water-filled  voids  within  the  sediment  fabric  (Rhoads  and 
Young.  1971). 

Micro  topography 

Bioturbation  produces  microtopography  of  the  .sea  floor.  Heezen  and 
Hollister  (1971)  have  suggested  that  the  visually  dominant  Ix-nthic  animals 
on  the  deep-sea  floor,  the  holothurians.  produce  “...  h'atures  mon'  wide¬ 
spread  and  more  visibly  evident  than  tho.se  produced  by  any  other  animal 
on  ('arth’’  by  their  extensive  sediment  reworking  activities.  Mounds,  depres¬ 
sions,  trails  and  burrows,  i.e.,  “Lebensspuren",  leading  to  pre.served  trace 
fossils  upon  lithification  have  been  subjected  to  intensive  study  for  decades 
(Schah^r,  1972;  Frey,  1975).  In  contrast,  the  importance  of  biogenicailly 
l)roduced  microtopography  in  generating  increased  bottom  roughness,  theia'- 
l)y  contriluiting  to  increa.sed  turbulence  at  the  sea  floor  and  increased  resus- 
pension  of  reworked  particles,  has  not  been  fully  apprc'ciated  (Rhoads.  lt)7  l  i. 

Bioturbation  activity  decrt'asc's  rapidly  below  tbe  uppc'r  10  cm  layer  of 
most  unconsolidated  marine  sediments  (Rhoads.  197  1).  alth.ough  la'rtain 
benthic  animals  have  b(>en  known  to  burrow  .several  mc'ters  deep  (.Myers, 

1979).  Di'pths  of  reworking  activity  in  deei)-s<‘a  sediments,  as  measured  by 
the  distribution  of  certain  radionuclides,  range  from  10  to  60  cm  (tluinasso 
and  Schink.  1975).  Since  liurrows  ;u-e  readily  recognizable  "Lebcmsspuri'n” 
they  are  commonly  recorded  from  deep-sea  cores  ( Chambi'rlain,  1975; 

Berger  el  al.,  1979).  F.lectron  microscopic  (>xaminations  of  epoxy-resin- 
impri'gnateti  Ihin-sections  of  continental-slo|)e  sinliments  havi'  dimionstrated 
clusti'ring  of  randomly  oriented  (larticles  and  vitv  largi'  voids  attributable  to 
extensive  burrfiwing  activity  (Bowles,  1968b).  Fvidence  of  quasifluid  defor¬ 
mation,  i.e..  false-bodied  thixotropy,  of  sediments  by  benthic  animals  is 
preserved  in  the  fossil  record  (Rhoads,  1970),  contrary  to  the  lontention  of 
Klliott  ( 1965).  False-bodied  thixotropy,  which  is  a  fluid  projierty  of  a  sedi¬ 
ment,  is  protluced  by  liiogenic  modifications  leading  to  decreases  in  com¬ 
paction  and  increases  in  porosity,  as  previously  discussed.  This  lyin'  of  tnogenic 
deformation  of  sediment  structure,  although  not  as  readily  recognizable  as 
that  resulting  from  burrowing  and  tube-building,  is  widespread  in  silty  mud 
sediments  and  may  be  responsible,  in  part,  for  the  open  cardhousi'  or  honey- 
coni  I)  sediment  fabric  reported  by  Bowles  ( 196 Ka). 
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PREDICTED  EFFECTS  OF  BIOTURBATION  ON  ACOUSTIC  PROPERTIES  OF 
SEDIMENTS 

It  has  been  shown  that  acoustic  properties  of  sediments  are  empirically 
predicted  as  well  as  theoretically  determined  by  physical  properties  of  sedi¬ 
ments.  It  has  also  been  shown  that  most,  if  not  all,  physical  properties  of 
sediments  which  are  important  to  geoacoustic  models  (e.g.,  porosity,  grain 
size,  density,  fabric,  rigidity  or  shear  modulus,  and  incompressibility  or  bulk 
modulus)  are  regulated  by  or  at  least  altered  by  bioturbation  by  benthic- 
animals.  We  therefore  hypothesize  that  bioturbation  by  benthic  animals 
alters  the  acoustic  properties  of  unconsolidated  marine  sediments. 

In  order  to  {iredict  the  magnitude  of  effects  of  bioturbation  on  marine 
si'diments,  we  calculated  the  affects  of  Nucula  annulata.  a  protobranch  bi¬ 
valve  (  Fig.l ),  on  acoustic  and  elastic  properties  of  sediments.  N.  annulata  is 
a  common  deposit-feeding  bivalve  found  at  station  “R”  in  Buzzards  Bay, 
Massachusetts  (Hampson,  1971).  The  sediment  at  station  “R”  is  primarily 
silt-  and  clay-sized  particles  (78—91%,  mean  grain  size  6  b)  as  determined  by 
standard  geological  grain-size  analysis  (Sanders,  1960).  The  intensive  sediment 
ri'working  by  N.  annulata  creates  an  uncompacted  granular  sediment  surface 
(  Fig. 2)  consisting  of  sand-sized  fecal  pellets  of  low  bulk  density  ( Young, 
1971),  The  ))orosity  of  near-surface  reworked'  station  “R”  .sediment  was  as 
high  as  86%  compared  to  70%  porosity  for  unreworked  station  “R"  sediment 
(  Rhoads  and  Young,  1971). 

The  velocity  of  compressional  waves  in  bioturbated  and  nonbioturbated 
sediments  was  calculated  from  Hamilton’s  (1974)  predictor  equation  for 
continental  terrace  sediments; 

V;  =  2455.9  -  21. 716(r;)  +  0.126(r?)-  (8) 

wliere  is  compressional  wave  velocity  (m/sec)  and  rj  is  porosity  (%).  The 
calculated  compressional  wave  velocities  were  then  corrected  for  the  20  m 
water  depth,  20°C  bottom  temperature  and  30° salinity  at  .station  “R" 
using  Wilson’s  ( 1962)  equations  for  the  speed  of  sound  in  sea  water 
(Hamilton,  1971).  Bioturbation  by  N.  annulata  increased  porosity,  which 
dei-reased  compressional  wave  velocity  from  1547  m/sec  to  1514  m/.sec 
(Table  I). 

The  attenuation  of  compre.ssional  waves  (a  =  dB/m)  in  bioturbated  and 
nonbioturbated  sediments  was  calculated  from  the  regression  equation: 

Q  =  0.7602  -  0.01487(rH  t  0.000078(tj  )Y  (9) 

where  rj  is  porosity  (%)  and  f  is  frequency  (8.5  Khz)  (Hamilton,  1972).  Com- 
t)ressional  waves  were  attenuated  less  in  bioturbated  sediment  than  in  sedi¬ 
ment  unaffected  by  bioturbation  (Table  I). 

Bottom  loss  (BL)  and  Rayleigh  reflection  coefficients  of  comjiressional 
waves  at  normal  incidence  to  the  sediment— water  interface  can  be  calculated 


■'  Th<'  term  '‘reworked"  is  used  synonymously  here  with  "hioturliuled". 
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Fig.l .  Sucula  annulata,  a  common  deposit  feeding  bivalve  found  at  Station  “R'’  in 
Buzzards  Bay,  Massachusetts. 


from  porosity -  impedance  relationships  given  by  Hamilton  (1970)  or  com- 
pressional  wave  velocity  times  density  formulae.  The  grajihical  and  comjnita- 
tional  values  were  the  same  for  eaeh  sediment  tyjie.  The  impt'dance  of  Itio- 
turbated  (1.89  gjem}  sec- 10®)  sediment  was  lower  than  nonliioturbated  (2.23 
g/cm*  sec  10®)  sediment.  The  Rayleigh  reflection  coefficient  calculated  from 
impedance  values  of  bioturbated  sediment  and  sea  water  was  0.10,  resulting  in 
a  20  dB  bottom  loss  for  compre.ssional  wave  at  normal  incidence  to  the  sedi¬ 
ment-water  interface.  The  sediment  which  was  not  bioturbated  had  a  calcu¬ 
lated  Rayleigh  reflection  coefficient  of  0.20  and  a  11  dB  bottom  loss. 
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A 


B 

Fig. 2.  A.  In  situ  photograph  of  reworked  sediment  at  Station  “R”  Buzzards  Bay, 
Massachusetts  (from  Young,  1971 ).  B.  Photograph  of  Buzzards  Bay  sediments  reworked 
by  Sucida  annulata  in  the  laboratory  (from  Young,  1971). 


The  incompressibility  (K)  or  bulk  modulus  was  estimated  using  Hamilton  s 
(1971a)  graphical  relationship  between  bulk  modulus  and  porosity  for  con¬ 
tinental  terrace  sediments.  The  bioturbated  sediment  had  a  lower  bulk 
modulus  and  was  therefore  more  easily  compre.ssed  than  sediment  which  was 
not  bioturbated  (Table  I).  The  rigidity  or  shear  modulus  (ju)  was  calculated 
from  eq.2.  The  bioturbated  .sediment  had  a  lower  calculated  rigidity  than  the 
nonbioturbated  sediment  (Table  I). 

The  shear  wave  velocity  ( V,)  was  calculated  from  eq.3  from  Hamilton 
( 1971a).  The  velocity  of  shear  waves  was  325  m/sec  in  nonbioturbated  sedi¬ 
ment  and  179  m/sec  in  bioturbated  sediment.  Although  not  calculated,  the 
attenuation  of  shear  waves  in  higher-porosity  bioturbated  sediment  would  be 
lower  than  the  attenuation  in  lower-porosity  nonbioturbated  sediment 
(Hamilton,  1976c). 
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TABLE  I 


Calculated  relationships  of  Station  “R”  (Buzzards  Bay,  Massachusetts)  sediments  bio- 
turbated  by  Nucitla  annutata  and  not  bioturbated  (porosity,  temperature  and  water  depth 
from  R  iioads  and  Young,  1970) 


Nonbioturhated 

Bioturbated 

Porosity 

im 

86% 

Temperature 

2WC 

20  C 

Water  depth 

20  m 

20  m 

Density  of  sediment  grains 

2.65  g/cmd 

2.65  g/cm* 

Density  of  water 

1.021  g/cm’ 

1.021  g/cm* 

Density  of  sediment 

1.51  g/cm^ 

1.25  g/cm’ 

Velocity  of  compressional  wave 

1547  m  'sec 

1514  m/sec 

Attenuation  at  3.5  khz 

0.35  dB/m 

0.20  dB/m 

Impedance  of  sediment 

2.33  g/'-m’  sec  •  10* 

1 .89  g/cm’  sec  -10' 

Rayleigh  reflection  coefficient 

0.20 

0.10 

Bottom  loss 

14  dB 

20  dB 

Bulk  modulus 

3.4  •  10'°  dynes/cm^ 

2.8  •  10“'  dynes/cm’ 

Shear  modulus 

0.16  •  10‘°  dynes/cm^ 

0.04  •  10'°  dynes/cm 

Velocity  of  shear  wave 

325  m/sec 

179  m/sec 

In  summary,  we  calculated  that  bioturbation  by  N.  annulata  dt'creased 
velocity  and  attenuation  of  compressional  and  shear  waves  in  sediments 
while  increasing  bottom  loss  at  the  sediment— water  interface.  We  then 
addressed  the  question:  “Do  the  calculated  changes  in  acoustic  and  elastic 
properties  of  sediment  bioturbated  by  N.  annulata  agree  with  the  theoretical 
relationships  between  physical,  elastic  and  acoustic  properties  of  sediment 
presented  by  Hamilton  (1971a)?’’ 

N.  annulata,  by  packaging  silt-  and  clay-sized  particles  into  sand-sized 
particles,  in  the  form  of  rod-like  fecal  pellets,  created  a  poorly  sorted  sedi¬ 
ment  with  a  lower  number  of  grain-to-grain  contacts  luid  a  more  open  fabric. 
This  firocess  resulted  in  lower  bulk  density  and  higher  |)orosity  of  sediment. 
The  sediment  would  be  expected  to  have  reduced  cohesion  because  of  the 
open  fabric  and  reduced  grain-to-grain  contact.  It  follows  from  Hamilton 
( 1971a)  that  reduced  cohesion  would  result  in  reduced  shear  modulus 
(rigidity)  of  bioturbated  sediment.  Rhoads  (197-1).  using  {H'netrometers  to 
measure  sediment  hardness,  confirmed  that  sediment  reworked  by  N.  annulata 
had  lower  shear  strength,  thus  lower  cohesion,  than  sediment  not  bioturbated 
by  N.  annulata. 

Although  bioturbation  can  change  characteristics  of  both  pore  water  and 
mineral  grains  (see  previous  section),  the  greatest  effects  of  bioturbation  on 
sediment  bulk  modulus  are  the  effects  on  frame  bulk  modulus  and  porosity. 
We  know  of  no  direct  or  indirect  measure  of  bulk  modulus  made  on  sediment 
bioturbated  by  N.  annulata.  It  is  logical  that  a  higher-porosity,  open-fabric, 
bioturbated  sediment  would  be  more  easily  compressed  than  a  lower-porosity 
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sediment  with  a  greater  numlier  of  grain-to-grain  contacts.  'I'he  poorly  sorted 
bioturbated  sediment  would  also  have  higher  permeability  and  decreased 
incompressibility. 

The  calculated  values  of  rigidity  and  incompressibility  and  the  measured 
value  of  bulk  density  agree  with  what  is  known  about  the  effects  of  bio- 
turbation  by  N.  annulata  on  sediments.  Since  the  acoustic  properties  of 
sediments  can  be  used  to  ciilculate  the  elastic  properties  of  sediments,  the 
pri'dicted  relationships  b('tween  bioturbation  and  the  elastic  and  acoustical 
properties  of  sediments  can  be  tested  by  measuring  the  velocity  of  comiires- 
sional  waves  in  bioturbated  and  unbioturbated  .sediments  in  laboratory 
e.speriments.  We  will  report  the  results  of  the.se  experimeiits  in  a  later  paper. 
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